ASYMPTOTICS OF UNIFORMLY RANDOM LOZENGE TILINGS 
OF POLYGONS. GAUSSIAN FREE FIELD 



LEONID PETROV 

Abstract. We study large-scale height fluctuations of random stepped sur- 
faces corresponding to uniformly random lozenge tilings of polygons on the 
triangular lattice. For a class of polygons (which allows arbitrarily large num- 
ber of sides), we show that these fluctuations are asymptotically governed by 
a Gaussian free (massless) field. This complements the similar result obtained 
by Kenyon IKen08l about tilings of regions without frozen facets of the limit 
shape. 

In our asymptotic analysis we use the explicit double contour integral for- 
mula for the determinantal correlation kernel of the model obtained previously 
in |Petl2l . 



1. Introduction and main result 

We begin with a description of the model and formulation of necessary previous 
results which motivate the main result of the present paper. The latter is stated in 
jjl~7| below. 

1.1. Model of uniformly random tilings. Consider a polygon drawn on the 
regular triangular lattice as shown on Fig. [T] 




FIGURE 1. Polygon on the triangular lattice. 

In the present paper we study the model of uniformly random tilings of such 
polygons by lozenges (= rhombi) of three types: 

Z7<>^ 

An example of such a tiling is presented on Fig. [2j left. Equivalent formulations of 
the model include: 

• (dimer interpretation) Lozenge tilings of a polygon are in a bijective correspon- 
dence with dimer coverings (= perfect matchings) on the part of the honeycomb 
graph located inside the polygon, see Fig. [2] right. 
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Figure 2. Tiling of a polygon on the regular triangular lattice 
(left) and its dimer interpretation (right). 

• (stepped surfaces) One can view each tiling such as on Fig. [2j left, as a 2- 
dimensional projection of a stepped surface, i.e., of a continuous 3-dimensional 
surface glued out of 1 x 1 x 1 boxes with sides parallel to three coordinate lines 
in space. 

The model of uniformly random lozenge tilings of polygons has received a sig- 
nificant attention over the past years }CKP01| . |KO07j . [Ken08] . See also [Ken09] 
for a detailed exposition of the subject and more references. 

1.2. Afflne transform and the class of polygons. For technical convenience, 
we perform a simple affine transform of lozenges which were present on Fig. [2] (see 
Fig. [3]). After this transform, tilings of polygons will look like the one on Fig. [7] 

Figure 3. Affine transform of lozenges. 

below. Polygons which are tiled will thus be drawn on the standard square grid, 
with all sides parallel either to one of the coordinate axes, or the vector (—1,1). 
We will denote the horizontal and the vertical integer coordinates on the square 
grid by x and n, respectively. 
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Figure 4. A polygon from the class we consider. In this example 
k = 4, and the polygon has 3fc = 12 sides. 

We will restrict ourselves to polygons of a special kind as shown on Fig. [4] 
Every polygon P we consider can be parametrized by two integers N = 1,2,..., 
and k = 2, 3, .. . (the polygon has 3fc sides), and by 2k (proper) half-integers 

Ax < B 1 < A 2 < B 2 < ... < A k < B k , A h Bi e Z' := Z + ±, 
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subject to the condition J2i=i(Bi — Ai) = N which ensures that there is at least 
one lozenge tiling of P. The bottom side of P lies on the horizontal axis n = 0, 
and all the k — 1 top sides (the ith such side has endpoints (B,, N) and (Ai + i, N), 
i = 1, . . . , k — 1) lie on one and the same line n — N . 

1.3. Height function of a tiling. We will view our lozenge tilings (as on Fig. [7] 

below) as projections of 3-dimensional stepped surfaces onto the [x, n) plane. The 
surface itself can thus be interpreted as a graph of a function h[x, n) which is called 
the height function of the tiling. To be concrete, let us stick to the convention that 
lozenges of type ^ correspond to horizontal planes, i.e., planes where the height 
function is constant. We also require that h(x,n) is zero near the lower left corner 
of the polygon. See Fig. [5] for an example and §3.1| below for a precise definition. 
See also |Ken09[ §2.8] for more discussion. 




Figure 5. Values of the height function h(x,n) on each horizontal plateau. 



The main object of the present paper is the height function corresponding to 
the uniformly random lozenge tiling of a polygon P. We assume that P belongs 
to the class of polygons described in j jl.2| We denote this random height function 
by hp(x,n). 

1.4. Limit shape. We consider large N asymptotics of random tilings as all di- 
mensions of the polygon P = P(iV) grow. That is, let the parameters Ai(N),Bi(N) 
of P(N) behave as 

Ai(N) = [diN] + const, B. t (N) = [b t N] + const (i = l,...,k). (1.1) 

Here a\ < b\ < . . . < dfc < are new continuous parameters which satisfy 
Si=i(^» — a i) = 1- The constants above are bounded uniformly in N and are 
neede d to ensu re that Ai(N),Bi(N) E U and £f =1 (-Bi(iV) - Ai(Nj) = N. 

In |CKP01] it was shown that in this N — > oo regime, the rescaled random 
stepped surface concentrates around a non-random limit shape which can be ob- 
tained as a unique solution to a suitable variational problem (see also [CLP98] . 
|DMB97j . and |Des98j ). This solution was described in |KO07j by means of the 
complex Burgers equation. For polygons we consider in the present paper, the 
limit shape is an algebraic surface. 

More precisely, the limit shape result means that the height function ftp(jv) obeys 
the following law of large numbers (with almost sure convergence) : 

h P (N)([xN},[ V N}) 

— N >h(x,V), N^oo, (1.2) 
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where (x,Tj) are the new global continuous coordinates, and h^x, 7 ]) is the function 
whose graph is the limit shape. The new coordinates (x, rj) are assumed to belong 
to the limiting polygon V which is paramet rized by {a^, in the same way as 

it was for P(N) and {A^N), B^N)}^ in fjl^J see Fig. [6] The new polygon V is 
located inside the strip < 77 < 1. 



«4 b/\ 




Figure 6. The limiting polygon V on the {Xi r l) plane and the 
frozen boundary curve. 



A feature of the model we deal with is that the limit shape develops frozen facets 
where the function h{x, ij) is linear. In other words, frozen facets correspond to 
zones inside V where lozenges of only one type are asymptotically present. Along 
with the frozen facets, there is a connected open liquid region T> <ZV . For (x, rj) £ 
T>, the limiting height function h^x^v) is curved: asymptotically inside the liquid 
region one sees a random mixture of all types of lozenges. E.g., see |CLP981 Fig. 2], 
BG09, Fig. 5], |KO07[ Fig. 1] for illustrations of uniformly random tilings with 
small mesh where the limit shape and the frozen boundary are clearly seen. 

Remark 1.1. There are more precise results in this direction. Namely, the as- 
ymptotic local distribution of lozenges around a given global position (x, rj) £ T> is 
governed by an ergodic translation invariant Gibbs measure on tilings of the whole 
plane. Such a measure is unique up to fixed proportions of lozenges of all types 
|She05| . and these proportions depend on the slope of the limit shape at the given 
point (x, rf). For polygons in the class described in j |1.2| this result was established 
in [Pet 12] , For k = 2 (when the polygon is a hexagon) this was obtained earlier 
in BKM M07j , |Gor08j . For tilings of regions when the limit shape has no frozen 
parts, the same local behavior was shown in [KcnO§]. See also |OR03| and [KOS06 
for more detail on the limiting translation invariant ergodic Gibbs measures. 

1.5. Complex structure on the limit shape surface. Complex coordinate on 
the limit shape surface was introduced in |KO07j (see also [Kcn08 ), and by a 
different technique in |Petl2j . From the results of [Pet 12] it follows (see j |4.2| below) 
that there exists a diffeomorphism w = w(x, rj) from the liquid region T> to the 
upper half plane i:= {zeC: 3z > 0}. (In |Petl2j the complex coordinate w was 
denoted by w c .) The function w(x,r?) is algebraic, it satisfies the following degree 
k equation: 

(w-x)!]* = (w-X+l-r?)TT fc >-M, (1-3) 
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and a version of the differential complex Burgers equation (see also |KO07| ): 

W (X,V) ~ X t dw(x,y) = dw{x,y) 
l — i] dx drj 



(1.4) 



The complex coordinate w(x, t]) can be used, in particular, to describe the local 
asymptotics of random tilings mentioned in Remark |1.1| (see |Petl2[ §2.3-2.4] for 
more detail). The complex structure w(x,?y) on the liquid region T> is employed in 
our description of asymptotics of fluctuations of the height function, see §1.7| below. 



1.6. Gaussian free field. Before we proceed to describing our results, let us first 
briefly discuss the object which governs the asymptotics of fluctuations of the height 
function, namely, the Gaussian free field. This subsection is adapted from |She07j . 
see that survey for a detailed and systematic discussion. 

The Gaussian free (massless) field GFF on the upper half plane H is a probability 
Gaussian measure supported on a suitable class of generalized functions on HI (and 
not on ordinary functions). In particular, the value GFF(z) at a point z £ H does 
not make sense. 

The distribution of GFF may be understood as follows. For any sequence of 
compactly supported smooth test functions {<j> r }^Li, the pairings {GFF(0 r )}^ 1 
form a sequence of mean zero Gaussian random variables with covariances 

E(GFF(fo)GFF(>,))= / \dz\ 2 {Vcj> k {z), V^(z)) 

\dzi\ 2 \dz-2,\ 2 (j>k{zi)4>i{z2) G(zi 1 z 2 ). 



Here (•, •) in the first integral is the usual inner product, and 

z,weH, (1.5) 



G{z,w) := In 
In 

is the Green function for the Laplace operator on the upper half plane EI with 
Dirichlet boundary conditions. 

Even though the value of GFF at a point cannot be defined, one can still think 
that the expectations of products of values of GFF at pairwise distinct points 
Z\ , . . . , z s are well-defined and are given by 

{s/2 
Ens(M«-i).M*)) seven ; 
0, s odd, 

with sum over all fixed point free involutions (= pairings) a on {1, . . . , s}. Indeed, 
for a finite number of test functions, 

E(GFF(0!)...GFF(&))= f E(GFF(z x ) . . . GFF(z 5 )) TT Idz^M^)- (1-6) 



The moments (1.6) uniquely determine the Gaussian free field. 
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1.7. Results. Now we are in a position to describe the main results of the present 
paper. We are interested in asymptotics of fluctuations 



H N (X,V) ■= h F{N) (\xN], [r,N]) -Eh P{N) ([ X N], [ V N]) 
of the height function (of a uniformly random tiling) around its mean. 



(1.7) 



Theorem 1.2 (Moment convergence of fluctuations to GFF). For pairwise distinct 
points (xi, T)i), . . . , (xs, 77 s ) inside the liquid region T>, as we scale the polygon P(N) 
as in the following convergence of moments holds: 



Jim ^ 2 E (H N ( X i,m) ■ ■ ■ H N (Xs, Vs)) 

= E ( GFF(w( X i, m)) ■ ■ ■ GFF(w( Xs , %))) 



(1.8) 



s/2 



^2~[[G (w(Xa(2i-i),Vv(2i-i)), w (Xa(2i),Va(2i))) s even; 

a i—l 

0, s odd, 

where the sum is taken over all fixed point free involutions a on {1, ... ,s}. 

Theorem 1.3 (Central limit theorem for fluctuations of the height function). The 

random function t/t{Hn{x, v) 011 "D weakly converges as N — > oo to the w-pullback 
of the Gaussian free field GFF on EL 



Theorem 1.3 means that for any smooth compactly supported test function </> on 
T>, we have the weak convergence as N — > oo: 

Hx,v)H N (x,v)dxdv -> / &(X,V) G ^( w (x,v)) d xdv (1-9) 
Jv Jv 

0(w~ 1 (z))J(z)GFF(z)|dz| 2 , 

where J(z) is the Jacobian of the change of variables z — > {XjV) by w _1 (which 
in fact can be explicitly calculated using ( 1.3 1-( 1.4|). Theorem |1.3| follows from 
Theorem |1.2| plus an additional bound on moments of fluctuations of the height 



function at infinitesimally close points (see \ 5A). 

It is worth noting that Gaussian free field fluctuations in random tiling models 
were also obtained (along with |Ken08j ) in |BF08| . |Duillj . and |Kuall| . 

1.8. Strategy of the proof and organization of the paper. Our proof is based 
on an explicit formula for the determinantal correlation kernel of uniformly random 
tilings which was established in |Petl2| . We recall these results in fj2] In f|3]we write 
the multipoint fluctuations E (Hn(xi, T)i) ■ ■ ■ Hn(Xs, Vs)) of the height function in 
terms of that correlation kernel. This allows us to establish Theorem 11.21 and then 
Theorem 1.3 in fusing certain fine asymptotic properties of the correlation kernel 
which are obtained in ^j4j 

Our argument generally follows the approach of BF08^ (especially see Section 5 
in that paper) which in turn was partly inspired by |Ken08| . We also use some ideas 
from [Duill] . However, our correlation kernel has a more complicated structure 
than those of [BF08] and |Duillj : the critical points of the action ( |4.2| which 
are solutions of (1.3) cannot be determined explicitly. Thus, in our fj4] in order to 
investigate asymptotics of the kernel we must employ certain new considerations 
QjjjOpH). 
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2. DETERMINANTAL STRUCTURE OF RANDOM TILINGS 



In this section we recall the formula of |Petl2| for the determinantal correlation 
kernel of our model of uniformly random tilings. Then we extend that kernel and 
describe the joint distribution of all three types of lozenges in our random tiling. 



Except for [2.4 this section is essentially taken from |Petl2j . 



2.1. Interlacing particle arrays. Let P be a polygon of our class in the (x, n) 
plane (see [1.2) contained inside the horizontal strip < n < N. We pass from 
tilings of P to interlacing particle arrays as follows. We first trivially extend any 
tiling of P to a tiling of the whole strip < n < N with N small triangles added on 
top (see Fig. [7]). Then we place a particle in the center of every lozenge of type ^. 




Figure 7. Tiling of a polygon and the corresponding interlacing 
particle array. 

Thus, we get a particle array X := {x™ : m = 1,...,N; j — 1, . . . , m} e 
2W(w+i)/2 w ith precisely m particles at each mth horizontal level, m = 0, 1, . . . , N. 
Because we started from a tiling of P, these particles must satisfy the interlacing 
constraints: 

xf +1 < x™- 1 < x? (2.1) 

(for all j's and m's for which these inequalities can be written out). Moreover, the 
particles in the top row of the array are fixed: 

{x% < . . . < xf } = {A x + \ < A 1 + § < . . . < B 1 ~ § < B l - 1 < (2.2) 
< A 2 + \ < ... < B 2 -\ < ... < A k + \ < ... < B k - \} 



(here {Ai, Bi}^ =1 are the parameters of P, see [ 1.2 1. Clearly, lozenge tilings of P 
and such interlacing arrays X with fixed top row (2.2) are in a bijective correspon- 
dence. For a connection of these arrays with Gelfand-Tsetlin schemes (an object 
related to branching of representations of unitary groups) see, e.g., |Petl21 §3]. 
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2.2. Determinantal correlation kernel. We see from §2.1| that the uniform mea- 
sure on the set of all tilings of the polygon P is the same as the uniform measure 



on the space of interlacing integer arrays X = {x™} with fixed top row (2.2 1. We 
denote both measures by Pp. Viewing X as a particle configuration, we can also 
think of the measure Pp as of a point process on Z x {1,2,..., N}. 

Definition 2.1. Let {x\, ni), . . . , (x s , n s ) be pairwise distinct positions, Xi € Z, 
1 < rii < N . The correlation functions of the point process Pp are defined as 

p s (xi , n\ ; . . . ; x s , n s ) :— Pp (there is a particle of the random configuration {x™} 

at position (xi, n{) for every i = 1, . . . , s). 

It is well-known that the measure Pp on interlacing particle arrays is determi- 
nantal (for instance, this fact can be deduced from the Kasteleyn theory, see j ]2.3| 
below). That is, there exists a function K(x,n;y,m) (the correlation kernel), such 
that 

p s {xi,ni\ . . .;x s ,n s ) = det[K(xi,ni;Xj,nj)] s ijj=1 (2.3) 

for any s and any collection of pairwise distinct positions (xi, ni), . . . , (x s , n s ). 
About determinantal point processes in general see the surveys [SosOOj . |HKP V06] . 
[Borllj . 

In [Pet 12] the following explicit formula for the correlation kernel K of random 



interlacing arrays X with fixed top row (2.2) was obtained 



Theorem 2.2 f |Pet!2p . For 1 < m < N, 1 < n 2 < N - 1, and x x ,x % G Z, the 
correlation kernel of the point process Pp has the forr^ 

(xi - x 2 + l) ni _„ 3 _i (iV-rii)! 



K(xi,ni;x 2 ,n 2 ) = -ln 2 <n 1 lx 2 <x 1 - 



(m-n 2 -l)! (JV-na-1)! 



1 / d f dw (z-x 2 + l) N ^ a - 1 1 A (Ai + 1 W ) Bi . Ai 

I ill) TiKr , i 11) 



(27ri) 2 / J (w- Xi) N _ ni+ i w- z Al + | _ z) B ._ A . ' 

£(2:2) c(oo) 

(2.4) 

The contours in z and w are positively (counter-clockwise) oriented and do not 
intersect. The contour £(x 2 ) in z encircles the integer points x 2 ,x 2 + l, . . . , — \ 
and only them (i.e., does not contain x 2 — l,x 2 — 2,... and B\~ + ^, Bk + §, . . ■). The 
contour c(oo) in w contains £(x 2 ) and all the points X\,X\ — 1, . . . ,X\ — (N — n\). 

The above explicit formula for the kernel K is our main tool in the present paper. 

2.3. Inverse Kasteleyn matrix. Here let us recall the connection |Petl21 §6] 
between the above kernel K (Theorem |2.2[ ) and the Kasteleyn matrix of the honey- 
comb graph Gp inside our polygon P (Fig. [2j right) . We will use it to write down 



below. 



the joint distribution of three types of lozenges ^, □, and W in f 2.4 

The honeycomb graph Gp is bipartite; its vertices correspond to two types of 
(triangle) faces in the dual triangular lattice: 



^Here and below denotes the indicator of a set, and (y) m '■= y(y + 1) • ■ ■ (y + m — 1), 

m = 1, 2, . . . (with (y)o := 1) is the Pochhammer symbol. 
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We will encode each such triangle by the position (x, n) of the mid-point of its 
horizontal side. The Kasteleyn matrix of the graph Gp is its adjacency matrix 
with rows and columns parametrized by white and black triangles, respectively 
(e.g., sec [Kcn09 ). Inside the polygon, this matrix looks as 

'1, if (y,m) = (x,n); 

,1, if (y,rn) = (x.n — 1): ,„ 
Kast(^(a;,n);W(y,m)) = { . ' ' (2-5) 

1, if (y, m) = (x + 1, n — 1); 

^0, otherwise 

(see Fig. [8|. For N(.x, n) on the boundary of the graph Gp, the ^(x, n)-th row of 



Figure 8. Edges of three directions in the graph Gp encoded by 
pairs of triangles. 

Kast will contain less than three ones, and the same for the k(y,m)-th column. 

It is known that the determinant det[Kast(^(a;, n); L.(y, to)], where ~^i(x,n) and 
L.(y,m) run over all possible white and black triangles in Gp, is equal to the total 
number of lozenge tilings of the polygon P. As |Ken09[ Cor. 3] suggests, Kast -1 
can serve as a correlation kernel for the uniform measure on tilings of P. A more 
precise statement looks as follows: 

Theorem 2.3 f |Petl2j ). The inverse Kasteleyn matrix and the correlation kernel 
K of Theorem \2.t\ are related as follows (for all possible values of (x, n) and (y, to) ): 

Kast _1 (k(j/,m);^l(a:,n)) = (-l) y - x+m - n K(x, n; y, to). 

2.4. Extension of K and joint distribution of three types of lozenges. 

Here we compute probabilities that a random tiling has lozenges of prescribed 
types at prescribed positions (i.e., the joint distribution of all types of lozenges). 
These probabilities are given by determinants similar to the correlation functions 



of interlacing particle arrays (Definition 2.1|, but with an extended kernel 



Let, by agreement, the position of every lozenge be encoded by the position of 
its white triangle (i.e., by the position of the circle dot on Fig. [sb. We introduce 
the following extended kernel (here and below 9j € {'si, □, ^.} are types of lozenges): 

K(xi,ni;x 2 ,n 2 ) 1 if 62 = n; 

—K(xx,ni;X2,ri2 — 1), if 6* 2 = (2.6) 



K s (x 1 , rai, 0i ; x 2 , n 2 , 9 2) ■= < 



K(xi,m;x2 + 1,U2 — 1), if 2 = 

Proposition 2.4. For any collection of lozenges of types 81,... ,8 S G {^, □, ^} at 

(pairwise distinct) positions (xi,m), i = 1, . . . , s, we have 

Pp ( There is a lozenge of type r at (x r , n r ) for all r = 1, . . . , s) 

= det [Kg (xi , n; , 0,- ; Xj , rij , 8j )] s itj=1 . 
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Proof. This is a direct consequence of the above Theorem 2.3 and also of [Ken09( 
Cor. 3]. Namely, 

• there is a lozenge of type ^ at (x, n) iff the dimer covering contains the edge 
(k(x,n);N(x,n)); 

• there is a lozenge of type □ at (x, n) iff the dimer covering contains the edge 
(k(x,n- l);^(x,n)); 

• there is a lozenge of type ^ at (x, n) iff the dimer covering contains the 
edge (k(x + l,n — 1); ^\l(x, n)); 

Then, using [Kcn09, Cor. 3], we can write the probability 

P P (There is a lozenge of type r at (x r , n r ) for all r = 1, . . . , s) 

as a determinant of a suitable matrix with entries Kast -1 . In this matrix there will 
be three types of rows (recall that they are indexed by black triangles) corresponding 
to different types of lozenges as above. Then, using the relation between Kast -1 

□ 



and the kernel K (Theorem 2.3), we conclude the proof 



A similar property for a different tiling model (of an infinite region) was obtained 
in [BF081 Thm. 5.2]. See also |BGR10l §7.2]. 

3. Height function and its multipoint fluctuations 

In this section we discuss the concept of a height function of a tiling. In |3.3| 
for our model of uniformly random tilings we express the multipoint moments of 
fluctuations of the height function through the correlation kernel K of Theorem|2.2| 



3.1. Definition of the height function. Let P be a polygon from our class 



({ 1.2 1. Fix a tiling of P. It is possible to define the height function of this tiling 
which at every position (x, n) € P is equal to 



h(x, n) 



l{there is a lozenge of type 'si or n at (x, m)}. (3.1) 



Clearly, this implies that the height function is constant on each horizontal plateau 
consisting of lozenges of type see Fig. [5] 

With every tiling one can associate three families of nonintersecting lattice paths 
(e.g., see [Pet 121 §2.5]). Nonintersecting paths in one of these families shown on 




Figure 9. "Level lines" of the height function (left) and the cor- 
responding parts of the frozen boundary (right). 

Fig. [9] (left) can serve as "level lines" of the height function. Namely, /i(x, n) at 
a given point is equal to the number of these nonintersecting paths lying between 
(x,?i) and the line n = 0. See also Fig. [7] where the tiling is extended so that 
formula (3.1 1 and the interpretation with the "level lines" make full sense. 
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Remark 3.1. The two other families of nonintersecting paths |Petl2[ §2.5] give 
two other possible ways to define the height function. See also |Ken09[ §2.8] for a 
more detailed discussion. 

3.2. Paths to the boundary. Now let P(iV) be a sequence of polygons scaled as 
explained in §1.4| We would like to study asymptotics of multipoint fluctuations of 
the random height function hptm (corresponding to the uniformly random tiling of 
P(JV)) around global positions • • ■ , (xs, %) € V. Here V is the liquid region 

(Fig. [6| . Namely, we are interested in the asymptotic behavior of the expectations 



entering (1.8 1 



E(H N ( X i,m)---H N (xs,Vs)), (3.2) 



where £/jv(Xj! %')' s are the fluctuations of the height function defined in (1.7| 



To compute values of the height function entering (3.2), one could use formula 



(3.1). But then several indicators corresponding to the same point (x,m) in (3.1) 
will enter the resulting expression. This will lead to certain technical complications 
which we can easily avoid. Namely, let us choose s piecewise-linear simple paths 
with horizontal and vertical pieces such that each ith path connects ([xiN,rjiN]) 
and the lower left part of the frozen boundary as on Fig. [10] We choose this part of 
dT> because on the corresponding facet the height function will be asymptotically 
equal to zero, cf. Fig. [5j By agreement, we assume that while crossing the frozen 
boundary, each path goes in vertical direction and proceeds vertically down until 
it hits the boundary of the polygon. By taking N larger if necessary, we require 
that these s paths do not intersect)^] We also require the number of piecewise-linear 
segments in each path to be bounded (uniformly in N). 




Figure 10. Paths to the frozen boundary along which we calcu- 
late the height function at ([xiN], [rjiN]), i — 1, . . . , s. 

Such paths can be constructed using the diffeomorphism w: T> - 



(see ( 1.5 



and (4.2 below) which maps the frozen boundary to the real line. It is not hard to 
show that in H the desired (continuous) nonintersecting paths exist. Then in T> we 
can approximate the images of these paths under w _1 by piecewise-linear paths. 
Since the number of points {xi, (and paths) is finite, we can also make sure that 
the number of segments in these paths is bounded uniformly in N. 

Remark 3.2. We could also use paths ending at any part of the frozen boundary 
on Fig. [9] (right) because on the corresponding facets our height function asymp- 
totically becomes constant, and we subtract this constant in the definition of Hn 



2 Except for the case when some of the points (xiiVi) coincide; then we still do not allow 
intersections away from the starting points ([XiiV], [fJjiV]). 
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(1.7). However, we use only paths as on Fig. 10 to simplify the notation in S4.6 



below. 

Definition 3.3. For integers x < x' and n, denote 

x 1 

^x,x'{n) := ljthere is a lozenge of type 'si at (y, n)}. (3.3) 

y— cc+1 

Also, for integers n < n' and y, set 

n 

V n , n '(y) := E l{there is a lozenge of type 'si or □ at (y,m)}. (3.4) 

m— n+l 

Above we have explained how each individual fluctuation Hisiixii ^7i) , * = 1 , - - • , s, 
can be written as a finite linear combination (with coefficients ±1) of expressions 
of the form H XtX i(ri) — EH X)a /(n) and V n ,ri{y) ~ EV n , n /(y)n Moreover, if all the 
points (xiiVi) are distinct, then each indicator entering any of the sums (3.3) and 
(3.4) corresponding to r)i), i = 1, . . . , s, appears only once because our paths to 



the boundary do not intersect. 

3.3. Expectation of a product of horizontal and vertical sums. Let all the 



points {Xh r ]i)i * = 1, ■ ■ ■ ,s, be distinct. From the discussion of S3.2 it follows that 



our expectation (3.2) can be expressed as a linear combination (with coefficients 

(3.5) 



±1) of terms of the form 

Em (H^jtnO-EH^nj)) f[ ( V «^(^) ~ EV «,,»J (xj)) 

M=l j = r+l 

such that the following horizontal and vertical segments 

{(y, Hi): y = Xi + l,..., x'j}, i = 1, . . . , r; 
{(xj,m): m = rij + l,.. . ,n'j}, j = r + 1,.. . ,s. 

do not intersect. Here and below we assume that x^ < x' t , rij < nf- for all 



Proposition 3.4 (cf. [BF08, Lemma 5.3]). With the above notation and assump- 
tions, expression \375V can be written in the following form: 



E •■• E E •■■ E ^t 

yi=X!+l y,=x r +l m r+1 =n r+1 + l m s =n s + l 



A.2,1 A2.2 



(3.6) 



The matrix blocks are given by 

Ai,i = [(1 - 5ij)if(yi,ni;yj,nj)]i,i=i,...,r 

Ai,2 = [K(yi,rii;xj + l,rrij - l)]i=i,..., r; j=r+i,...,s 

-4-2,1 = [-K(xi,mi;yj,nj)] i=r+1> _ >s .j = t t ... ir 

A 2 ,2 = [-(1 - 5ij)K(xi,mi;xj + l,m,j - l)]i,j= r +i,...,s- 



(3.7) 



Note that the height function vanishes at the end of each path on Fig. 10 
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Proof. As was observed in |Ken08l Proof of Thm. 7.2] , the subtraction of the means 
in (3.5) leads to vanishing of the diagonal matrix elements in Ai i and A 2j 2 in (3.7). 
Thus, it suffices to consider E ( Y\[=i H Xi ,a:< i n i) E[j=r+i ^nj , < ( x j )) ■ We write every 
H Xi!X /(nj) and V n ^ n '.(xj) as the corresponding sum (3.3H3.4). After taking the 
expectation, we get an s-fold sum as in (3.6) with terms 

Pp(AT) (There is a lozenge of type *sl at (yi, rii), i = 1, . . . , r; 

and of type 'si or □ at (xj^rrij), j = r + 1, . . . , s) 

= I"p(Af) (There is a lozenge of type 'si at (yi, n^), i = 1, . . . , r; (3.8) 

and there is no lozenge of type ^ at (xj , irij ) , j = r + 1, . . . , s) . 

The latter probability can be expressed as an s x s determinant which has the 
structure 



det 



— l,...,r; j ,j' — 



Here is the kernel which is obtained from Kg ( 2.6 ) via a particle-hole involution 
(e.g., see BOO00, Appendix A. 3]) at positions {(xj,rrij)} s j =r+1 . That kernel Kg- 
looks as follows (see also Proposition 2.4): 

^(yi,™^;yi',»V,^) = K{y i ,n i ;y i >,n il ); 
K^yi^rii^Xj, ,m r ,■%.) = ^(j/^ni,^;^/,™^,^) 

= K(yi,ni\Xj> + l,rrij> - 1); 
K^Xj.m^^yi^rii,^) = -K e (xj, ra 3 ■, y v , n v , ^) = -K{x^m^y v ,n v ); 



Kg (Xj , TOj , IEj' , TTij 



= djji — K(xj,rrij] xy + 1, my — 1). 



Then, setting the diagonal matrix elements to zero, we obtain matrix blocks (3.7). 
This concludes the proof. □ 

Remark 3.5. Representation of Proposition |3.4| is not valid if some of the points 



(XjiVj) coincide because then we cannot write all the probabilities (3.8) as s x s 
determinants. In this case, in the asymptotic analysis of multipoint fluctuations 
(3.2) we employ Lemma 5.6 below. 



4. ASYMPTOTICS OF THE KERNEL 

In this section we investigate asymptotic properties of the correlation kernel of 
Theorem |2.2| in various regimes. 



4.1. Action S(w; x, i])- The polygon P(-/V) is assumed to be scaled as in £ 1.4 We 
will be interested in asymptotics of the kernel K(x\, m; X2, ni) when the two points 
{x\,n\) and (x2,n 2 ) behave as 



N 



-> Xj, 



N 



1,2, 



(4.1) 



where (xiiVi) arL d {X2tV2) are two (not necessarily distinct) global positions inside 
the limiting polygon V (see \ 1.4 and especially Fig. [6|. 
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Definition 4.1 ( !Petl2[ §7.2]). Define the action by 

S(w; x, rj) := (w - x) ln(w - x) - {w - x + 1 - v) m ( w ~ X + 1 - v) ( 4 - 2 ) 

k 

+ (1 — 77) ln(l — rj) + — w) ln(6j — w) — (a,; — u>) ln(<2j — w) . 

i=l 

Unless otherwise stated, we assume that that the branches of all logarithms have 
cuts looking in negative direction along the real line. Note that the real part 
3tS(w; x, v) is well-defined and continuous for all w G C. 



Denote also 



E(w;x,r)) 



(w - x)(w -x + i-y) 

1 — 7] 



and 



Ej(w) := E(w; %, Sj(w) := S(w; U, ^f), 



j = l,2. 



(4.3) 



Proposition 4.2. In the regime (4.1), the kernel K(xi,rii;x 2 ,n 2 ) of Theorem 2.2 
has the following asymptotics: 



K(xi,nx;x 2 , n 2 ) = -l„ 2<ni (l + O(^)) 



1 
2~7ri 



e X p{N(S 1 (z)-S 2 (z))} 



dz 



+ (l + 0(^)) 



(27Ti) 2 



e'(X2-) ^i(z)S 2 (z) 

1 exp{JV(5i(ii;)-5 2 («))} 



dz <z> du> 

C( X 2-) JcCoo) «>-z A /Hi(ty)5 2 (z) 



(4.4) 



-ffere z m both single and double integrals runs over a counter-clockwise contour 
which crosses the real line just to the left ofx2, and also to the right ofbk ~ ^jf- (see 
(1.1)). The w contour is counter-clockwise, contains £(x2— ) (without intersecting 
it), and is sufficiently large. 

When {xiiVi) — (X2>V2)> this essentially coincides with |Petl2i Prop. 7.2]. 



Proof. Let us adapt the double contour integral in formula (2.4) to the asymptotic 
regime (4.1 1 by scaling the variables as z — z/N, w = w/N (and then renaming 
back to z, w): 

v, \ i i (x! - x 2 + l) ni -n 2 -i . N(N-ni)\ 
K{Xi,ni;X2,n 2 ) = -l„ 2<ni l X2 < Xl ; — h — — x 



(2m) 



^2 



dz 



dw 



(ni - n 2 - 1)! 

(Nz - X 2 + l)jv-na-l 1 



C(X2 



c(oo) 



(Nw - a;i)jv-m+i w - z ^ [Ai + \ - Nz) Bi 



{N-n 2 -l)V 
(Ai + i - Nw) Bt -A i 



Here z and w run over the corresponding scaled contours, they can be c hosen 
independent of TVF] These contours coincide with the ones in the claim (4.4). 



4 We may drag the z contour slightly to the left of X2 because the integrand has zeroes in z 



which allow that; and also drag it to the right of because the integrand in ( |2.4| does not have 
z poles to the right of — | . 
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Expressing all the Pochhammer symbols in the integrand above through the 
Gamma function and applying the Stirling approximation, we may write for non- 
real z,w (see |Petl2[ §7.2] for more detail): 

1 N(N-m)\ (Nz-X2 + l)N- n2 -i -rT {A i + ^-Nw) Bi -A i 
w-z(N-n 2 -l)\ (Nw-x 1 ) N - ni+1 y (A i + \-Nz) Bi -A i 

= (l + O(jr)) / exp {N(S 1 (w) - S 2 (z)) } . 



As for the additional summand, using Lemma 6.2 in Petl2], we write 

1 -, ( x l — %2 + l)»ii-n 2 -l 



L Ti2 <ni ^X^<X\ 



(ni -n 2 - 1)! 



(4.5) 



= _ 1 {N~ ni )\ x J_f (Z-X 2 + l)jV-n 2 -l ^ 

" 2<ni (_/v - n 2 - 1)! 27ri/ c(x2 ) (z-x^n-^+i 

Then, scaling the z variable as above for the double integral (z — z/N), and using 
the fact that 

N(N- ni )l (Nz-x 2 + l) N - ni -i , n( 1 ^ exp{jV(^i(z)-^ 2 (z))} 
(JV-na-1)! (Nz- Xl ) N _ ni+ i ( + [n>) ^E 1 (z)E 2 (z) 

we conclude the proof. □ 



4.2. Critical points of the action. Proposition |4.2| suggests to use the saddle 
point (steepest descent) approach (e.g., see |Oko02[ §3]) to investigate the asymp- 
totics of the correlation kernel K(xi, ni; x 2 , n 2 ). The first step is to understand 
critical points of the action, i.e., points where S'(w; %, rj) := ^S(w,Xy r l) = 0- Let 
us recall the results about critical points obtained in |Petl2| : 

1. Depending on the global position (x, rj) inside the limiting polygon V, there are 
either or 1 critical points of the action in the (open) upper half plane H. 

2. Points (x,v) £ V for which there exists a nonreal critical point (denote it by 
w (x,v)) constitute the (open) liquid region V c V where asymptotically one 
sees all three types of lozenges (see Remark |1.1| ) . 

3. As a function of the global position (x, rj) G T>, w(x, rj) satisfies the algebraic 



equation ( IL3J and a form of the complex Burgers equation (1.4). 
4. When (x, rj) £ V approaches the frozen boundary curve dV (which separates 
the liquid region from frozen facets), the critical point w(x,ry) € EI merges with 
its complex conjugate w(x, rj). In addition, points (x,v) € dT> that are cusps 
(= turning points), or points where &D is tangent to a side of the polygon (see 
Fig. [6]) correspond to certain more special types of merging of the critical points 
w(x, rj) and w(%, rj), which wc do not need to address in our analysis]^] 

Thus, for all (x,v) € 9T>, the action S(w;x7 r l) nas ( at least) double criti- 
cal point w(x, rj) £ K which can be taken as a real parameter on the frozen 
boundary curve. The map w _1 : M — > dV is one-to-one and rational (see [Pet 121 
Prop. 2.6]). 

Proposition 4.3. The map w: T> — > H, (XiV) ^ w{x,vj), * s a diffeomorphism. 



^See the explanation in the proof of Lemma 4.15 that one can choose paths on Fig. |lo| away 



from such more special points. 
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Proof. Finding the image of a point z£i under the inverse map w 1 amounts to 
solving the equation JO} 



(*-x)TT. Az-a l ) = (z-x + l-ri)\^. (z-h) (4.6) 

for % and 77. Since 2 € EI is complex and (x, 77) must be real, this is actually a pair 
of real equations. Let us first rewrite (|4.6|) as 



*-« P-IL,^ +&+"- 1 )IL, 



Since the imaginary part of nJ=i Ib^T * s nonzero for z 6 H (Lemma 4.4 below), 
one can always solve the equation 

for % + 77 — 1 and thus find a solution (x, 77) which belongs to T> because we have a 
bijection on the boundary dT>. This implies that the map w : T> — > H is bijective. 
The map w is differentiable, and, moreover, the partial derivatives w x and 



cannot both be zero inside V because of the complex Burgers equation (1.4 1. One 



can also check that the inverse map is differentiable. This concludes the proof. □ 

Lemma 4.4. Let a\ < b\ < . . . 

of the limiting polygon V . Then 



Lemma 4.4. Let a\ < b\ < . . . < < bk, 53j=i(^i — a i) ~ 1> ^ e ^ e parameters 



It, 



z 



z — Oj 

■—hi 



^0, ze 



Proof. Observe that the argument of * * is the angle under which the segment 

[oj, hi] is seen from the point z. Thus, the argument of the whole product Y\.%=\ t-a* 
must be strictly between and tt, and so the imaginary part of that product cannot 
vanish. □ 

4.3. Moving the contours. Our aim in this subsection is to explain how we 



deform the contours in the double integral in (4.4) to employ the saddle point 
analysis. 

Let us assume that (not necessarily distinct) limiting global positions (xiiVi) 



and (x2,f]2) m (4.1) belong to the liquid region V C V. Denote the corresponding 
critical points of the action by Wj := w (%-, ^-j £ H, j = 1, 2. 

The behavior of Si. 2 around W1.2 is quadratic because these critical points are 
simple (see also the proof of Proposition 4.11[ ). Thus, there are four curves starting 



from each point w 12 along which the imaginary part S£(Si %{w)) is constant (see 



Fig. 11). As the new w contour we choose the counter-clockwise closed contour 
passing through Wi composed of two curves with $s(Si(w) — Si(wi)) = on which 
$t(Si(w) — Si(wi)) < for w 7^ wi,wi (Fig. [IT] left). The new counter-clockwise z 
contour must pass through W2 and look like the one on Fig. [TlJ right, so on it we 
will have ^(S^z) — ^(wa)) > for z^W2,W2. 



Proposition 4.5. The z and w contours in the double integral in (4-4) can always 



be deformed to become the new contours described above (indicated on Fig. 11) 



ASYMPTOTICS OF RANDOM LOZENGE TILINGS. GAUSSIAN FREE FIELD 



17 









/z 




t).4 
02 


( ( | 


\\ 




OD 
- ()2 








in 


A 



Figure 11. Critical points wi (left) and w 2 (right). Along the 
bold curves one has $s(Sj(-) — Sj(vjj)) = 0, j = 1,2. Shaded are 
regions where 5i(Sj(-) — Sj(wj)) > 0. The new w and z contours 
are indicated on the left and on the right, respectively. 



This results in the following asymptotics of the kernel K in the regime G77 



K(x 1 ,n 1 ;x 2 ,n 2 ) = (l + O(^)) 



1 

2tti 



dz 



exp{N(S 1 (z)-S 2 (z))} 



(4.7) 



(2rri) 5 



V5!(z)H 2 (z) 
dzdu; exp {iV(5iH - } 
{2} w - 2 1 /Si(w)S 2 (z) 



where in the double integral {z} and {w} are the new deformed contours. 

The single integral may or may not be present; this depends on whether the new 
contours intersect, and also on the inequality between n\ and n 2 . All these cases 
can be unified by choosing an appropriate contour £ S ingio ; see §4-4\ 



Proof. Let us fix any (%, rj) £ T>, and set w := w(x, Tj) and S(z) :— S(z;x,v)- 
As our first step, we aim to justify that the picture of shaded regions where 
$t(S(z) — S(w)) > looks exactly as on Fig. 11 and also describe the points where 
the four contours {z: ^sS(z) = 35(w)} intersect the real line. 

Because diS(z;x,n) ~ 77 In 1 2; | as \z\ — > 00, and < n < 1 inside V, far away on 
Fig. 11 we see a shaded region, i.e., where $t(S(z) — 5(w)) > 0. 

Since S(z) is holomorphic everywhere in H, along each of the four contours 
{z: ^SS(z) — SJS'(w)} (the thick curves on Fig. [Tlj) the sign of $l(S(z) — S(w)) must 



be constant. This implies that each thick curve on Fig. [TTJfrom w to w must be 
completely inside a shaded or not shaded region. 

Now let us look at the function ?s(S(z) — SS'(w)) for z £ K + ie for fixed small 
e > 0. Observe that 

+ ie) ln(i + ie)) = eln \t + ie| + t arg(t + ie) ~ tt • (t)_ := rr ■ tl t < , 

where t € K, as e -» 0+. Thus, 



1 c> 



9f5(t + ie) ~ (t - x)- ~ (t - X + 1 - V)- + - =1 K ft i - *)- - ^ - *)-]' 



see Fig.pj Clearly, ^^(i + ie) ~ l te[x+ ,_ 1)X] - £)* =1 hel^-,^]- Looking at the 
slopes ofthe graph on Fig. fl2| we see that the four contours {z: SsS(z) = SJS'(w)} 
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FIGURE 12. Graph of ±$S(S(t + ie)-S(w)), t e M, for small e > 0. 
The segments [aj,bj] (red) and [x + V^ 1;X] (black) are displayed. 



can intersect the real line in at most three points)^] Because of the relation between 
these contours and the shaded regions on Fig. [TT] explained above, there are exactly 
three such points of intersection: 

• t+ € [x + V - 1,%], where 9fS(t+) = 3S(w) and ftS(t+) > ftS(w); 

• < t~ , both belonging to the union of the segments [aj,bj], where ^sS(t7 _) = 
35(w) and K5(tr.) < »5(w). 



Moreover, from Fig. 12 we see that t, < t + < t r . The fourth contour {z : QS(z) 
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SS(w)} (with $iS(z) > %iS{w)) runs to infinity, see Fig. 

Now as a second step, we explain how we can move the z and w contours in 
the double contour integral in (4.4) to get (4.7). Looking at the poles in z and w 
in the original integrand in (2.4), we see that: 

• We can drag the points of intersection of the z contour with R (without picking 



any residues) everywhere except in regions where the slope of the graph on Fig. 12 
is strictly negative. 

• The same for the w contour: we cannot drag it through regions where the slope 
of the graph on Fig. 12 is strictly positive. 

The old z and w contours are described in Proposition |4.2[ together with what 
was said above, we see that these z and w contours can always be deformed in a 
desired way. The new w contour will intersect the real line at points i ; ~~ r (xi, f)i)', 
the new z contour — at i + (x2,7?2)- Because the integrand in (2.4) is regular in z 
at z — oo, we can let the z contour pass through infinity. 

In the course of this deformation no residues coming from poles on K will be 
picked. However, if the new z and w contours intersect, the residue at w — z will 
be picked from the w integral, and then this residue will be integrated in z over an 
appropriate arc. Together with the single integral already present in (4.4), this will 
lead to appearance of the single integral in (4.7). We will describe and investigate 
it in 34.41 below. □ 



In fact, the case when there are infinitely many such points (i.e., when a horizontal part of 
the graph on Fig. |12| is lying at the horizontal coordinate line) corresponds to (XiV) belonging to 
the frozen boundary, when the action has a real double critical point (see [4.2 k 
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4.4. Estimating the single integral. The goal of this subsection is to asymp- 



totically estimate the single integral in (4.7 1. A priori from the proof of Proposi 
tion 



4.5 we see that the integral over £ s ingic may look as follows (we omit the factor 

exp{N(S 1 (z)-S 2 (z))} dz y 



1 + O(jj)) and the integrand 



l 

27ri 



(a) If the new z and w contours of Proposition 4.5 do not intersect, then the 
integral has the form l„ 2 > ni §{ z }, where {z} is the full new z contour. 

(b) If the new contours intersect at points £ G EI and then the integral has the 



form — 1 



n 2 <n 



MO 



ln 2 >m j£ R (ni wnere £l(C) is the left part of the new z 



contour passed from £ to £, and Cr(C) is its right part passed from £ through 
oo to £. 

See Fig. [13] for the possible configurations of contours. 




Figure 13. Various possibilities for the new deformed z and it) contours. 

Remark 4.6. If the new z and w contours intersect more than once in H, then the 
contour £ s ingio would contain several parts. However, then we always can write an 
estimate of the form 

/ | . . . | dz < l n2 <m / \...\dz + l 1l2 > ni \...\dz 
•Mingic JZhiO J<£r(C) 

(dots mean the integrand), where C and C' ar e some points of intersection of the 
new contours. Below (in Lemmas 4.8 and 4.9 and Proposition 4.10[ ) we estimate 
the above two summands separately, so we may think that the case (b) covers all 
possibilities when the two contours intersect. 

First, we deal with the case (a): 

Lemma 4.7. If the new z and w contours do not intersect, then there is in fact 



no single integral in (4-1) 



alent to (see Proposition 4.2) 
1 



Proof. Recall that the single integral in (4.7) in that case is asymptotically equiv- 

l)A r -n 2 -l 



1 



n 2 >rii 



27ri 



iV(JV-7ii)! (Nz 



x 2 



7i2-l)! (Nz-nt + 1) 



-dz. 
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Here {z} is the full new z contour. This integral can be explicitly computed using 
Lemma 6.2 in |Petl2j (see also (4.5)), it is equal to 



1 



(ni 



n 2 



(x 1 -x 2 )\ 

This expression is nonzero only if x\ + n\ < x 2 + n 2 , otherwise the Pochhammer 
symbol vanishes. But observe that the three inequalities 



n 2 >ni, x 2 <xi, 



xi + ni < x 2 + n 2 



in the regime (4.1 ) imply that (for large TV) the segment [x2+?72 — 1, X2] is completely 
inside [xi + 771 — 1 , xi]. From the proof of Proposition 4.5 it follows that the new 
z contour crosses the real line at some point inside [x 2 + 7y 2 — 1, X2] (and hence 
inside [xi + 771 — 1, xi])i an d the new w contour passes through two real points at 
the opposite sides of [xi + t]i ~ 1; Xi]- Thus, we see that in this situation the new 
z and w contours must intersect. This concludes the proof. □ 

Now we will obtain certain estimates in the case (b). 

Let n 2 > ni, and £r(C) f or C € H be defined as above. We have the 



Lemma 4.8 

estimate 

1 



2?ri 



eMNiSAz) - S 2 (z))} 



<£«(£) V /S l( z )"2(z) 



dz 



<Cexp{N-$t(S 1 (0-S 2 (0)}- 



Here the constant C is uniform for (x\, Hi), (x 2 , n 2 ) in the regime (4-D with the 
condition n 2 > n\, for the limiting global positions (xi)?7l); (X2)^2) belonging to a 
compact region T> c d T> . 

Proof. Assume first that {x\,n{) ^ {x 2 ,n 2 ). For large \z\, we have the expansion 
F(z) := Si(z) - S 2 (z) = const + \nz 



2 N 2 



(1 



f)(l 



U2_\ 

N I 



(l-f)(l-f))+0(i). 



Observe that the function F(z) has no nonreal critical points. This implies that 
there is a curve in EI starting at the point ( e H along which 3F(z) = 3^(0 and 
$iF(z) < $tF(£) for z 7^ (. This curve can either extend to infinity, or cross the real 
line somewhere in the segment [X2 + V2 — LX2] (see Fig. 14). This can be seen by 
considering the function QF(t +ie) of t € M similarly to the proof of Proposition 4.5 
Note that for n 2 = ni, such curves will never to go to infinity (Fig. 14 right) 
the lower half plane the situation is symmetric. 

Since the integrand is regular at z = 00 (see the proof of Proposition 4.2), we 



In 



can always transform the contour £r(C) so that it will consist of curves described 
above (along which 5RF(z) < $IF(()). This implies the claim for (x\, n\) ^ (x 2 , n 2 ) 
because if {XjiVj) S T> c , then the factor 1/ y / Ei(z)'E, 2 (z) in the integrand is uni- 
formly bounded. For (xi,ni) — (x 2 ,n 2 ), the integral does not depend on N, and 
the claim is trivial. □ 



Lemma 4.9. Let n 2 < n\, and £l(C) for £ s 
1 / eMNiSxiz) - S 2 (z))} 



2ni 



Ct(C) v/^i(z)S 2 (z) 



dz 



be defined as above. We have 
<Cexp{N-^(S 1 (0-S 2 (0)}. 
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Figure 14. Level lines of ^RF(z) (top; warmer colors represent 
larger values) and QF(z) (bottom) in the case n 2 > ri\ (left) and 
n 2 = n\ (right). The red (left) segment is [xi + Vi ~ liXi]> an d 
the black (right) one is [x2 + m - 1,X2]- 



Here the constant C is uniform for (x\, n\), (aca, 71%) in the regime (4-1) with the 
condition n 2 < n\, for the limiting global positions {xi,Vi)i (X2, T l2) belonging to a 
compact region T> c <g T>. 



Proof. This is established in the same way as Lemma 4.8 Since (C) never extends 



to infinity, we can always transform it to get the desired estimate, see Fig. |15| □ 

From Lemmas |4.8| and |4.9| we derive a stronger estimate which we will use: 
Proposition 4.10. The single integral in \4- can be estimated as 



1 

2tti 



exp{N(S 1 (z)-S 2 (z))} 
A /S 1 (z)5 2 (z) 



dz 



< 



Cexp{N- ^(Q - S 2 (Q)} 
1 + R 



(4.8) 



'Single 

where R := \f (xi 

and (x 2 ,n 2 ) behaving as in (4-1), for the limiting global positions (xi,?7i), {X21V2) 
belonging to a compact region V c <s V. 



x 2 ) 2 + (^1 — n 2 ) 2 . The constant C is uniform for (xx,ni) 
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Proof. The passage from the estimates of Lemmas 4.8 and 4.9 to ( 4.8 1 is done simi- 



larly to [Duilll Lemma 6.3] and is based on a standard steepest descent argument. 
Let £ € HI be the point of intersection of the z and w contours of Proposition 4.5 



(see also Fig. 13 1 where the contour £ s ingie starts. (If the z and w contours do not 
intersect, the claim is trivial by Lemma 4.7 ) Assume that we have transformed 

so that on it we have 3F(z) = ^SF(() and ^StF(z) < 
62(2). If the new contour extends to 



4.8 



or 



4.9 



^single as in Lemma 

SRF(C) for z ^ C, where Fjz) = Si(z) 
infinity (Fig. 14 left), let us close is so that it will become bounded. Denote this 
new contour by Cg itlglc . 

Let us choose a smooth parametrization z = z(t) of the part of the curve £g ing i e 
inside H such that z(0) = £ and z(l) e M. We have \z'(t)\ < const (and by 
compactness of T> c this constant is uniform in (xj,rij)), so 



1 

27ri 



exp{iV(5 1 (z)-5 2 (z))} 



dz 



< const - / exp{7V • $t{F{z{t)))}dt. 
'a 



'Ctagi. ^B. 1 (z)E 2 (z) 
Let us assume that (x\, Tlx) ^ (a: 2 , n 2 ) otherwise the claim is again trivial. We have 
[x\ — X2, ni ~ n 2) = -R(cos 0, sin cf>), R > 0. 

Denote 

N-$t(F(z(t))) = R-G lj2 (t;<t>). 

The property that 3lF(z) < 3tF(£) on our contour allows to write the following 
estimate. Since V c is compact and ( depends continuously on (Xj,?7j), we can 
choose r > uniformly so that 

Gi, 2 (^)-G 1>2 (0;<£) <-At, 0<t<r, 

with a constant A not depending on £ or <f>. 
Then we have 

e RGi, 2 {t;4>) dt _ e flGi, 2 (O;0) f e fl(Gi, 2 (t;0)-Gi, 2 (O;0))^ 

I! JO 
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< e flGi. 2 (O;0) 



-ARt 



dt < —e 
AR 



RG lt2 (0;<t>) 



and 



e 



because on the contour {z(t) : r < t < 1} we have the inequalities ^tF(z) < 
$tF(z(r)) < 5RF(C) - Ar. This concludes the proof. □ 

4.5. Asymptotics of the kernel in the bulk. Our aim in this subsection is to 
write an exact asymptotic expansion of the kernel K(xi,ni;x 2 ,n 2 ) ( |4.7[ ) when the 
two points (xi,ni) and (x 2 ,n 2 ) are in the bulk of the system (i.e., they behave as 
in (4.1) and (XjiVj) € T>) and are sufficiently far from each other. The technique 
of getting such an expansion involves only "local" properties of the double contour 



integral formula ( |4.7[ ) for the kernel (in contrast to some considerations of f 4.4 >, 
and mainly follows the approach of [BF08] and [DuilT]. 



Proposition 4.11. Fix (sufficiently small) 5 > and a compact T> c <e T>. Uni- 
formly in (xj,nj) (j — 1,2) with € T> c , such that ||(xi, Tii) — (a; 2 ,n 2 )|| > 
N? +s , we have the following expansion: 



K(xi,ni;x 2 ,n 2 ) 
1 / 



,iV(Si(wi)-S 2 (w 2 )) 



2ttN\ ( Wl -w 2 ) A /Hi(wi)H 2 (w 2 )(-5i'(wi))5(5^(w 2 ))3 



3 Ar(Si(wi)-S 2 (w2)) 



(wi - w 2 ) v / S 1 (w 1 )H 2 (w 2 )(-S'i / (wi)) 2 (S 2 '(w 2 )) 2 

e iV(S 1 (w 1 )-S 2 (W 2 )) 

(wi - w 2 )^Si(wi)S 2 (w 2 )(-S , ('(wi))5(S , 2 '(w 2 ))2 

e Af(Si(w 1 )-S 2 (W 2 )) 

(Wi - w 2 )^S 1 (w 1 )S 2 (w 2 )(-S'i(w 1 ))2-(5 2 '(w 2 ))5 



(4.9) 



(l + 0(iV-*)). 



The branches of the square roots (iS 1 "^ • • )) 5 above are chosen in accordance with 
the directions of the w and z contours in the double integral in (^.1) at points W\, 
Wi 7 and w 2 , w 2 , respectively (see (4- ID in the proof). 



Proof. Observe that the contribution from the single integral in (4.7) given in 
Proposition |4 . 1 0] is asymptotically negligible in comparison with the desired expan- 
sion (4.9). Thus, it suffices to consider only the double contour integral in (4.7): 



I 2 {xi,ni;x 2 ,n 2 ) 



(27Ti) 2 



{z}J{w} 



dzdwexp{N(S 1 (w)~ S 2 (z))} 
- z v / E 1 (w)E 2 (z) 



(4.10) 



Recall that the w contour passes through the critical points Wi , Wi , and on it 
we have 5sSi(w) = 35i(wi), and $lSi(w) < 9?Si(wi) for w ^ wi,wi. Similarly 
for the z contour: it passes through w 2 ,w 2 , and on it ^sS 2 (z) = 3S' 2 (w 2 ), and 
5R5 2 (z) > 3?5 2 (w 2 ) for z ^ w 2 w 2 

The main contributions to (4.10) come from neighborhoods of the critical points; 
and parts of the contours which are sufficiently far from them give an exponentially 
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small contribution. Since there are four pairs of critical points, we get four sum- 
mands in (4.9). Let us consider only the case of (w!,w 2 ), for the other three pairs 
the situation is analogous. 

In small neighborhoods of Wi and W2 let us replace the (curved) w and z contours 
by the corresponding tangent lines. Introduce the local coordinates t, s G [— N 5 , N s ] 
as follows: 

t s 
w(t)=w 1 + ^= r , z(s) = w 2 + —= r . (4.11) 

Here the branches of the square roots (-SJ'(wi))3 and (5^'(w 2 ))3 are chosen so 
that when t (resp. s) increases, the point w(t) (resp. z{s)) passes along the tangent 
line to the w (resp. z) contour in the direction of that contour. 

In these new variables the behavior of the exponents in the double contour 
integral is 

lim N{S l {w{t)) - SiK)) = -t 2 /2, lim N(S 2 (z(s)) - S 2 (w 2 )) = s 2 /2. 

N— >oo JV— s-oo 

(4.12) 

The convergence here is uniform for t, s £ [— N 5 ,N S ] and also (by compactness of 
T> c and continuity) for our values of (xj, rij). Moreover, at the endpoints t, s = ±N S 
the expressions e Ar ( s i('"( t ))- s 'i( w i)) and e N(s 2 (z(s))-s 2 (w 2 )) are exponent i a ii y sma n. 
Parts of the contours which are even farther from the critical points Wi , Wi and 
w 2 ,w 2 thus give an exponentially negligible contribution. 

This implies that the double contour integral (4.10) picks the following contri- 
bution from the neighborhood of (w 1 ,w 2 ): 

1 1 f N ' f N " dsdt _ex P {N(S 1 (w(t))-S 2 (z( s )))} 



(27ri) 2 Ar(_^'( Wl ))l(^'( W2 ))l J_ N s J_ N s ' («,(*) - z(s)) A /Hi(w(i))H 2 (z(s)) ' 

(4.13) 

Let us now get rid of non-exponential terms in the integral above. The map 
w _1 : EI — > T> is a diffeomorphism, so there exists a constant A > such that 

|wi -w 2 | > ~ (X2,V2)\\ > AN-^ +s . 

The second derivatives S" 2 (wi j2 ) are nonzero inside T> and hence are bounded from 
below in T> c (recall that they vanish on the frozen boundary dT>). Thus, we may 
write 

,J n = + o{N-->)) 

w[t) — z{s) wi — w 2 

where the constant in 0(N~^) does not depend on S (it depends only on V c ). We 
can also replace Si(u>(t))S 2 (z(s)) by 5i(wi)S 2 (w 2 ), this will affect the asymptotics 
by a factor which is less significant than (1 + 0{N^^)). Thus, we may rewrite 
|4T3| ) as 

I e JV(Si(wi)-S a (w 2 )) 

" W? iV(-5^(wi)) 5 (S 2 '(w 2 ))* (wi - w 2 ) v /sT(w 1 )S 2 (w 2 ) X 

x (l + 0(iV-i)) / / e ^(Sx(«(t))-Si( Wl )-s 3 (*W)+&(w a )) dsdi> 



-JV S J-N 6 
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Taking N large and using the uniform asymptotics (4.12), we see that the above 
double integral becomes Gaussian and can be explicitly evaluated. This concludes 
the proof. □ 



Corollary 4.12. In the setting of Pr oposition J^.ll, we have the same expansion 
for K(xi, n\\ %i + 1, n 2 — 1) as in (4-9), but with an extra factor of ^E^m~ i n each 
term with w 2 , and with a factor of '^E^jjf ^ n each term with w 2 . 
Proof. This is obtained in the same way as Proposition |4. 1 3~| using the fact that 



See also Lemma 7.4 in [Pet 12] . 



)-l n (l-«a)+ln(z-f) 



□ 



We can also write an estimate of the double contour integral I 2 (4.101 when the 
points (xi,ni) and (x 2 ,n 2 ) are sufficiently close: 

Lemma 4.13. Fix (sufficiently small) S > and a compact T> c <s T> . Uniformly in 
(xj,nj) (j = 1,2) with (tv~)Tv) € D C) such that ||(xi,ni) — (^2 , ^2) 11 < N? +s , we 
have the following estimate: 

( J e N-U(S 1 {w 1 )-S2(v J2 )) 

\I 2 (x 1 ,n 1 ;X2,n 2 )\ < 



N 



Proof. We argue as in the above Proposition 4.11 but now we must estimate the 
term l/(w(t) — z(s)) in a different way. Since the points Wi and W2 are close, we can 
write S" (wi) = S 2 ' (w 2 )(l+0(A^5+< 5 )), where the constant in 0(N~^ +5 ) is uniform 
in our (xj,nj)'s and depends only on T> c . This implies that in l/(w(t) — z(s)) we 
can replace (— S"(wi))i with ±i(5 , 2 ' (w 2 ))5 , where the sign ± depends on the choice 
of square roots. Moreover, we have |wi — w 2 | = 0(A r_ 3+ (5 ) ; so we can write 



N 



w{t) - z{s) y/N( Wl - w 2 ) - (^(w 2 ))-»(« ± it) ' 

Then we proceed as in the proof of Proposition |4.11| and see that the resulting 
double integral has the following asymptotics coming from the neighborhood of 
(wx,w 2 ): 

1 pW(Sl(wi)-S 2 (W2)) 



(2tt)2 V^V(-5^(w 1 ))5(S^(w 2 ))VSi(wi)5 2 (w 2 ) 

N s r N s Jv(s 1 ( lU (t))-Si(w 1 )-S 2 ( 2 i( s ))+S 2 (w 2 )) 



-dsdt. 

i~n s J-n 6 VN(wi — w 2 ) — (S' 2 ' (w 2 ))~5 (s ± it) 

(For other three pairs of critical points, one can get the same estimate.) Depending 
on how close the points wi and w/2 in our regime, we sec that the above integral 
may have a singularity which is integrable, and (on the other hand) the expression 
\/N (vji — w 2 ) may go to infinity. This implies that we can always bound the integral 
by a constant, and thus we arrive at the desired estimate. □ 



Remark 4.14. In Proposition 4.11 we see that when the points (xj, rij) (j = 1, 2) 



are sufficiently far from each other, the main contribution to K(xi, ni;a; 2 ,n 2 ) (4.7) 
comes from the double contour integral. On the contrary, for sufficiently close 
points (xj,nj), the single integral in (4.7) is asymptotically more significant. On 
the extreme, when asymptotically the differences x\ — x 2 , ri\ — n 2 € Z stabilize, the 
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double integral in (4.7) vanishes in the limit, while the single integral gives rise to 



the incomplete beta kernel, see |Petl21 Thm. 2 and Prop. 7.9] 
4.6. Estimates of the kernel close to the edge and in the facet. We conclude 



this section with several estimates for the kernel K{x\, n\\ x 2 , n 2 ) (4.7) when one or 



both points (xj,rij) becomes close to the lower left edge of the liquid region V, or 
outside T> in the lower left facet, see Fig. [10] We mainly follow a similar treatment 
for a simpler kernel which was performed in [BFQEI §6]. 

Recall (?3.2 1 that we choose the paths for calculating the height function as 
on Fig. 10 ending in the lower left facet. Let 77 = ?7fb(x) be the equation for 



the corresponding lower left part of the frozen boundary. Thus, the liquid region 
(sufficiently close to that part of dT>) is determined by the inequality 77 — ?7ft>(x) > 0. 
We distinguish three regimes for a point (x, n) (in the pre- limit integer coordinates): 



(inside the liquid region) 
(close to the edge) 
(at the edge or in the facet) 



n-Nrjfoix/N) > N^; (4.14) 
TVs > n - iV?7f b (x/N) > cN^; (4.15) 
n-Nrjfbix/N) < cN* , (4.16) 



for some c > 0. 



Lemma 4.15. Assume that the points (xj,rij) (j — 1,2) behave as in (4--1), and 
one or both of them is close to the lower left edge as in \4.15 ). Also, let |wx — w 2 | be 



bounded away from zero uniformly in N . Then there exists c in (4-15) large enough 
so that we have 

C . e AfK(Si(w 1 )-S 2 (w 2 )) 

\K{x x ,m;x 2 ,n 2 )\< w =, (4.17) 

AV|SY(wi)S£(w 2 )| 

uniformly in N. 

Proof. Because |wi — W2I must be bounded from below, we see that the limiting 
global positions (XjiVj)) 3 = 1>2) must be distinct. Proposition |4. 10 (cf. Remark 
4.14) then implies that the single integral in (4.7) is asymptotically less significant 



than the desired estimate (4.17) for the kernel (note that at least one of the factors 
S" (wi), S'2 (w 2 ) goes to zero as TV — > 00, see also the proof of Lemma 5.2 ). Therefore, 
it suffices to derive (4.17) for the double contour integral in (4.7). 



As usual in the steepest descent approach, the main contribution to the double 



contour integral in (4.7) comes from the neighborhoods of the critical points. Thus, 
there we have w « wi, z ss W2 (plus three more possibilities with Wi^ replacing 
W1.2, but they give the same contribution to the desired bound). 



Let, by agreement, the paths on Fig. 10 be separated (in the limiting coordinates 



(x, v)) form the tangent points of the frozen boundary and the sides of the polygon. 
Clearly, such paths still can be chosen. Thus, we may think that the quantities 
Si(wi) and 2 2 (w 2 ) are uniformly bounded away from zero (see also |Petl21 Prop. 
2.7 and Fig. 14]). Thus, it remains to estimate the product of two single integrals 



e NS ^dwi e~ NS ^dz, 

{w} J{z} 



where the w and z contours are as in (4.7). We will derive the estimate of the form 

s-t iiVJRSx 2 (wi 2) 

—, for each of the single integrals (with "+" sign for the first integral, 

V JV l S l,2( w l,2)l 

and "— " for the second one), and this will give the desired claim. If, say, (xi,ni) 
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is not close to the edge, then the corresponding estimate can be obtained in the 
same way as in Proposition 4.11 So let us assume that (xi, ni) is close to the lower 
left edge, and estimate the w integral above; for the z integral the argument is the 
same. 

From [Petl2[ Prop 2.7] it follows that for ( £1,711 ) approaching the lower left edge 
of the liquid region, the corresponding critical point Wi approaches the real line to 
the left of the point a 1 from flllj ). Using |Petl2[ §2.3 and (2.10)] (cf. Remark [Ll) 
it can be shown that arg S"(wi) tends to — 7r/2. Let us introduce the local variable 
t around wi: 

w{t) = wi + e-' m/i t, -5<t< V2$s(wi). 

It can be readily checked (cf. |BF081 Lemma 6.8]) that replacing the w contour 
around Wi by the straight line {w(t)} will not affect the desired bound (provided 



that c in (4.15) is large enough). We then have 



o(t 4 



and 5R(— ^"(wi)^ 2 ) w — i|5"(wi)|i 2 . Since Wi is close to the real line, one can 
derive an equivalence of the form 5"'(wi) « 

l_j_j ^111 / N . 1 * I 



|Si(wi)| 

a(wi) ■ 



and so 



|^(wi)|t d 



6 v / 2S(wi) 

We see that for —6 < t < 0, the cubic term helps the convergence, and for < 
t < \Z23(wi), we can bound the cubic term by the quadratic term which will also 
ensure the convergence. We thus see that the integral of e iV (' 5 i(" , )- s i( w i)) around 
W! is equal to a constant times the integral of exp(— ~ \S"(wi)\t 2 ), which leads to 



the desired estimate. 



□ 



To describe further estimates, we need to introduce some notation. Let (x, n) 



be at the lower left edge or in the corresponding facet as in (4.16). We would like 
to mimic the critical point w (-^, ^) and the value of the action S(w(jj, ^); j|, j^) 
for such (x, n) as follows: 

w = w(f,£) :=v.(f A (f)) e i, 



S(i 



' AT' N 



k 

+ (1 - f ) In (1 - rm, (f )) + - w ) Hh ~w)-(a 



l 



!? + !-%>(#)) 

- w) ln(ai — w) . 



Denote by Wi^ and S\^(w) the corresponding quantities at (^1,2,^1,2) similarly to 
wi j2 and Si^iw) (see also (4.3)). Note that when the point (x,n) is at the edge, 
i.e., when n — Nrjo-, (x/N) 



and same for S. 



O(Ns), we have |w - w| = 0(Ni 

Lemma 4.16. Assume that the situation is as in Lemma \4-15\ but now the point 
(x2,n 2 ) is at the lower left edge or in the corresponding facet as in (^.16), while 



(a;i,m) is in the bulk or close to the edge (4-- 15). Let |wi — W2I be bounded 



away from zero uniformly in N . Then, also uniformly in N , we get the following 
estimate with C, C 2 > 0: 



\K(x 1 ,n 1 ;x 2 ,n 2 )\ < 



Cr e iV3ftSi(wi) 
y/N\S?(*i)\ 



-AfsRS 2 (vv 2 ) 



exp -C* 2 iVf k I)" f 



>} 
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Proof. Similarly to the beginning of the proof of Lemma |"4.15[ we see that it suffices 
to estimate the product of two single integrals. The w integral is bounded as in 
Lemma 4.15 (yielding the first factor in the claim). Thus, it remains to estimate 
the z integral §^ e~ NS2< - z ^dz. We will mainly follow the approach of |BF08[ §6.1] 
(which is in turn based on the technique first applied in BFS08, Prop. 15 and 17]). 
We provide a brief derivation omitting certain bounds which are done in a way 
similar to what is performed in [BF08, §6.1]. 

Let us first consider the following scaling of (x 2 ,n 2 ): 

X2 = IX2N}, n 2 = [Nift, (x 2 ) + uN%], 

where n£l, and \i is some coordinate such that the line {x- X = X2} intersects 
the lower left part of the frozen boundary as on Fig. |10| Let us expand 

-NS(z; f , ft ) „ -NS(z; X 2,Vfb (X2)) - uN^S v (z; X 2,Vfb (xa)). 
We deform the z contour in e~ NS ' 2 ^ dz so that it will pass through the real 
double critical point w 2 = W (X2 5 Vib (X2))- As it usually happens for Airy-type 
asymptotics, the main contribution to the integral comes from an N 3 -neighborhood 
of the double critical point. Let us introduce the local variable t: 

z = w 2 + tN~s, 

and continue the above expansion: 

- NS (w 2 + tN~ s ; X2 , Vfh (Xa)) - uN^> S v (w 2 + tN~ ^; X2 , Vfb (Xa)) 

rj -NS(w 2 ,X2,Vfb (x 2 )) - f\t 3 S / "(w 2 ;x2,Vib (xa)) 

- uNis ri (w 2 ;x2,Vfb (Xa)) - u*^(w 2 ;X2,»7fb (xa)) 

The terms - g t 3 S'" (w 2 ; %2 , ?7fb (X2)) - utS' v (w 2 ; Xa , »?fb (X2)) after the integration in 
the neighborhood of the double critical point w 2 give the Airy function (similarly to 
BF08, Lemma 6.1]). The terms containing N contribute to the factor e - N ^s 2 (w2) 
(after substituting u = N% ( ft - r) & (ft)))- 

For general (x 2 ,n 2 ) as in (4.16), the desired exponential estimate containing 
e const « £ or ^ e gj n gi e integral is obtained along the lines of Lemma 6.2 in [BF08J. 
This concludes the proof. □ 

Lemma 4.17. Assume that now both the points [xj, rij) (j — 1,2) are at the lower 
left edge or in the corresponding facet as in (4. 16). Let |wi — w 2 | be bounded away 
from zero uniformly in N . Then 

£ e AT(9^Si(wi)-?R,§2(w 2 )) 

\K(x 1 ,n 1 ;x 2l n 2 )\ < 5 x 



x expj-diVl ( m (ft) - ft) - C 2 7Vf (r ?fb (ft) - ft)} 
uniformly in N, where C,C\,C 2 > 0. 

Proof. This lemma is obtained similarly to the previous Lemma |4.16| but now we 
derive exponential estimates for both w and z integrals. □ 

5. Completing the proofs 



In this section wc finish the proof of Theorem 1.2 (§f 5.1-5.31, and then explain 
how it leads to Theorem 1.3 ({5A\. 
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5.1. Expanding determinants in s-fold sums (|3.6| ). Fix pairwise distinct points 
(Xii Vi)i • ■ ■ i (Xsj Vs) inside the liquid regio n T>. In 331 we showed that the expectation 
E (Hn(xi, Tji) . . . Hn(xs, V s)) of Theorem 1.2 can be expressed as a linear combi- 



nation of expressions like (3.6) 



E ■ 

y 1 =x 1 + l 



X , "r+l 

■ E E ■ 

y r =x r -\- 1 m r _)_i — n r -j_i + 1 



• E 

?7i s — n 5 +l 



det 



Ai, 2 

A 2 ,2 



(5.1) 



where the matrix blocks are given in ( |3.7| . Each such s-fold sum corresponds to a 
choice of one linear (horizontal or vertical) part on every ith path starting at the 
point {xi,Vi)i i = l,.-.,s (see Fig. [l0|. Throughout the section we assume that 
these paths on Fig. [10] along which we calculate the height function are separated 
from tangent points as explained in the proof of Lemma |4.15| 

Let us consider one s-fold sum as in ( 5. 1 ) . Expanding the above s x s determinant, 



we write it as the sum over permutations a G S(s) of terms each of which is sgnu 
times the product of matrix elements with indices (i, o~i), i = 1, . . . , s. Express a as 

"Aii A 1<2 ] 



a union of several disjoint cycles. Since the matrix 



-4, 



An 



has zero diagonal 



entries, all these cycles must have length > 2. In the next subsection we will show 
that the contribution of permutations containing cycles of length > 3 becomes 
negligible in the limit as N — > oo. 

5.2. Contribution of permutations with cycles of length > 3. Let the per- 
mutation o~ E £>(s) contain a cycle of length £ > 3. To shorten the notation, we 
assume that this cycle is 1 — >• 2 — >• . ..—>£—>!. In the expansion of the de- 



terminant in (5.1) we take the product of the kernels and do a horizontal (over 
Uj = Xj + 1, . . . , x'j) or vertical (over rrij = Hj + 1, . . . , n'-) summation. Let us 
collect terms corresponding to a fixed index i = 1, . . . , I. We will assume that the 
shifts i ± 1 are given mod I. There are four possible cases we consider: 

(V) The summation related to the index i is performed over a vertical segment: 
{(xi, rrii) : rrii = rii + 1, . . . , n'i} . It can happen that this vertical segment crosses 



the lower left frozen boundary (see Fig. 10). Thus, we split the summation into 



three parts according to ( 4.14 )-( 4.16 ). 



(V; edge or facet) Summation over Ii 



{rii + 1, ... ,7i^} n {mi-, mi < 
We need to consider 



Nrjfb (xi/N) + cNs}. Here r]^ is defined in [4.6 

K(ti-i,Ui-i;xi + l,rrn- l)K(xi,mi;ti + i,u i+ i) 



E 



(5.2) 



The minus sign is coming from the second factor (see (3.7)). Here and below the 
points (ti±i,Ui±i) (corresponding to indices Ui and a~ ) are equal to (yj.rij) or 
(xj + l,m,j — 1) for suitable j, see (3.7). 

Lemma 5.1. The contribution of the sum \5.!tfy over I\ goes to zero as N — > oo. 



Proof. By Lemmas 4.16 and 4.17 we can write 

E K(ti-x,Ui-i;Xi + 1,771,- - l)K(xi,mi\t i+1 ,u i+ i) 

E m . e/l ex p{- const • N ^ (r/fb ~ w )} X terms in (*i±l>W»±l) 



const 

< 



const . 

< j- x terms m (t i ± 1 ,u i ±i). 

TV 3 
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To get the first estimate above we employ considerations similar to Corollary |4.12[ 
this may change the bound only by a factor of a constant. The second estimate 
concludes the proof. □ 



(V; close to edge) Summation over I2 
A^fb (xi/N) < m t < Ni + Nr,^ (xi/N)}. 



{n, + !,...,<} n {m z : cN 3 + 



Lemma 5.2. The contribution of the sum (5.2), where rrii runs over I2 instead of 
I\, also goes to zero as N — > 00. 



Proof. We have 



E 

const 

< 



K(ti-i,Ui_i;Xi + l,m 
1 



E 

rrii£l2 



V ^ |SJ'(w(i))| 
Here we have used a shorthand notation w(i) 



l)if(af i) m i ;ti + i,«i + i) 
x terms in (i»±i, Ui±i). 



w (/v"' Tv )' anc ^ same f° r Si (see 



also (4.3 1). Also, as in the previous lemma we use argument similar to CorollarylT 



but this again may only change the constant in the bound. 

Let us bound |S'"(w(i))| from below. Observe that the third derivative S"'(w) is 
bounded away from zero for small 5u>, if %iw belongs to a bounded interval to the 
left of the point a\ <G K from (1.1 ). Indeed, under the map w _1 : H — > T>, such w's 



are close to the lower left edge and are separated from the tangent points. Thus, 
we can write |)Sf(w(i))| > const • 3w(i). 

Next, it can be seen that the imaginary part 3w(i) can be bounded from below by 
const • |w^(x, n)\ ■ |^ — ?7fb(ff ) + 0(N~®)\, where (x, rj) € V is some intermediate 
point closer to the edge than , r jf-)- Instead of 0(N~^) one could take any 
correction term which is asymptotically smaller than ^ — ?7fb(fr )• We have (see 
also [PetHl §7.6]): 

w-x 



E(w) 



E 



k 



r) — 1 + (w — x)( w — X + 1 — v)^( w ) 

where 

1 1 

= 1 V W — &i w — fl,; 

Using the formula for the action ( |4.2[ ), we see that 

1 



w 



Therefore, 



|S7(w(z))| > const 



S"(w(x, v);x> v) • ( w - x + 1 - v) 
Vfb (f) + O(N-h) 



(5.3) 



N 



S"(w(x,ri);x,v) 



> const • N 3 



N 



Vfb (f) + 0(N-°) 



(we estimate the second derivative in the denominator using (4.15)). We obtain 



const 



N ^ 



1 



|Sf(w(i))| 
This concludes the proof. 



< 



E 



const • N 



m t -N Vfh (j-) + 0(NB)\ 



< const • N~3 In AT. 



□ 
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It is not hard to see from Lemmas |5.1| and |5.2| that in the two remaining cases 
(V; bulk) and (H) (see below) we may assume that all the I variables we are 
summing over (which correspond to vertical and horizontal summations) always 



belong to the bulk of the system in the sense of (4.14) 



(V; bulk) Summation over I 3 :— {m+1, . . . , n^}n{m i : to, > As +Nr) fb (xj/iV)}. 
From what was said right above, we may as well take I3 = {rii + 1, . . . , n^}. We 
will investigate the asymptotics of ( |5.2| w here now runs over ^3 instead of I\. 

Similarly to the proof of Lemma 5.2 let us denote w(i) = w(^ , ^), w(i ± 
1) := w(^ i , ^r 1 ), and same for Si,Ei (see also (4.3)). Also, let /3 z {i) denote thi 
argument of the tangent vector to the z contour at the point w(i) as on Fig. 11 
and analogously for f3 w (i). It can be readily checked (using the global structure of 
the z and w contours (j|4.3[) that 



2wq (for some q £ Z). 

(5.4) 



&(*) = M%) + f , and 0„(i) + f3 w (i) + argS|'(w(i)) 



By the nature of our paths (on Fig. 10), the points (ti±i,Ui±i) are sufficiently 
far from (jXi,mi). Thus, we may use Proposition 4.11 and Corollary 4.12 and write 



(here and below S > is sufficiently small and fixed) 



E 



mie/3 



K(ti-i,Ui-i;Xi + l,m,i- l)K(x l ,m i - 1 t l+1 ,Ui + i) 



l + 0{N-2) 



2ttN 



E 



m,eza |£?(w(i))| 

JMi) 



1 w(i)-Xi/N 



w(i — 1) — w(i) w(i) — vj(i + 1) 3j(w(i)) 1 — rtii/N 

e^W e-'ft-W e -2JVi-ag i (w(«)) w( ,) _ 

w(i — 1) — w(i) w(i) — w(« + 1) |Sj(w(i))| 1 — rrii/N 

w(i — 1) — w(£) w(i) — w(i + 1) |S,(w(i))| 1 — rrii/N 

e -i/3*« e -ife(i) 1 w(l)-Xi/N 

w(i — 1) — w(j) w(i) — w(i + 1) S,(w(i)) 1 — rrii/N 



(5.5) 



x (terms in (tj±i, itj±i)) + O 



Here and below O denotes all additional terms (there are 12 of them in the above 
formula) which are obtained by replacing w(i — 1) and/or w(i + 1) by the corre- 
sponding complex conjugate points. 

Arguing analogously to Section 5.3 in BF08] (case (a/3)), one can show that 
the contribution of the oscillating terms above containing e ± 2Ari ' 3 ' s i( w ( 1 )) becomes 
negligible (of order 0(A~s +c )) in the limit. The remaining terms are smooth, and 
change over distances rrii ~ A. Therefore, up to an error of order 0(N~3), we 
can replace the summation over in ( |5.5[ ) by integration in the scaled variables. 
Namely, setting fi := m,i/N, r\ := (rij + 1)/A, r)' := n'jN, \ '-~ X i/N, we can 
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rewrite (5.5) as 



i + ow-f) P , 1 



2tt 



W0)l 



1 w(i) - x 



w(i — 1) — w(i) w(i) — w(i + 1) 5j(w(i)) 1 — /j 

e -i/3»(i) e -i/3„(i) ! w(i)-X 

w(? — 1) — w(i) w(i) — w(i + 1) Si(w(i)) 1 — /i 



(5.6) 



x (terms in (t i±ll u i±1 )) + Q j. 

The next step we perform is a change of variables. For the term with w(i), we set 
Cj + := w(i) = w(x,/i). The integration path r+ for C 4 . is from w(x,f?) to w(%, 77'), 
i.e., T+ is the image of the vertical line from 77 to 77' in T> under the map w : T> — > H. 



Form the results of |Petl2| . we have (see also (5.3)) 

d(j~ ^ 1 w(z) - x 

dp, S'! (w(i)) • Si(w(»)) 1-M ' 

Symmetrically, for the term with w(i), let £r := w(i) = w(x, /i), and the integration 
path r~ for C~ is conjugate to that for C z + - It can be readily verified (in particu- 
lar, using (5.4)) that the above integral (5.6| can be rewritten as the sum of two 
integrals: 



1 

2vri 



E 



1 



w(» - 1) - cr cr - w(i + 1) 



x terms in (t l± i,u»±i) + 



(5.7) 



There is an additional minus sign for IV because of a different phase —(/3 z (i)+l3 w (i)) 



in the second summand in (5.6). Thus, we have established the following fact: 



Proposition 5.3. The sum (5.2), where m>i runs over I 3 instead of I±, is (up to 

a factor of 1 + 0(N~z ), where S > is a fixed sufficiently small constant) equal to 
the sum of two integrals (5.1). 



(H) The summation related to the index i is performed over a horizontal seg- 
ment: {(y^m): yi =Xi + l,.. .jE'J. 

Proposition 5.4. The horizontal sum has the following asymptotics: 



E 



s , % l + 0(7V-s) 
K(t i ^ 1 ,Ui-i]y i ,ni)K(y i ,n i - 1 t i+1 ,u i+1 ) = — x 

yi=Xi + l ZTT1 



E 1 

Ej=± 



1 



1 



w(< - 1) - cr cr - w(» + 1) 



x terms in (U±i, Ui±i) + O 



where all the notation is as in the previous case, except that now the path of inte- 



gration connects w ^^ 1 , ji' 



and w(?f, j^), and T i is the conjugate ofTf. 



Proof. This is established in the same way as Proposition |5. 3 



□ 

After considering the four above cases, we conclude this subsection with the 
desired statement about permutations cr € <5(s) having cycles of length i > 3: 
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Proposition 5.5. Consider one s-fold sum (5.1) and expand the s x s determinant 
as a sum over permutations a £ ©(s). Then the contribution of those o~'s having 
cycles of length £ > 3 goes to zero as N — > oo . 

Proof. From Lemmas 5.1| 5.2 and Propositions 5.3 5.4 it follows that each cycle 
ji — > j2 ■ ■ ■ — ► ji — > ji in o~ asymptotically produces the following sum of £-fold 
integrals: 



1 



(27ri)' 



■ E 

Ei,...,e e =± 



ei... et 



i 



On the other hand, by |Ken08[ Lemma 7.3], we have 

s 

1 



all ^-cycles r£S(<) i=l Tl 



n 



0. 



This concludes the proof. 



£> 3. 



□ 



5.3. Contribution of fixed-point-free involutions. In §S 5.1 5.2 we have shown 
that if one expands the determinant in (5.1 1 as a sum over permutations a G ©( s )> 
then the contribution of permutations a which are not fixed-point-free involutions 
(— p airings) becomes negligible in the limit. Collecting all summands of the form 
(5.1 ) corresponding to the expectation E (i?jv(xi; Vi) • • • Hn(Xs, %)) (with pairwise 
distinct positions (xi, f]i), . . . , (xs, %)), we see that 



E (H N ( X i,Vi) ■ ■ ■ H N ( X s, Vs)) = (1 + 0(JV"5)) 

"(Xor(2i-l)l'7<r(2i-l)) 



(X<r(2i-l),»7 c r(2»-l)) 



S/2 

e n 

pairings a 6 &(s) i= 

w(X< r (2i).»7<r(2i)) 

77 ' 

(,t,2i-l — S>2iJ 



1 

1 



(x 



(2i).'7 CT (2i)J 



Note the additional minus sign coming from the signature of each transposition. 
The paths of integration in H from vj(xa(j)T r la(j)) to w(x CT (j)) Va(j)) are obtained 
by linearity (£3.2) and by symmetry of the contours in Propositions 
Each integral above can be explicitly evaluated: 



5.3 



5.4 



1 



(27Ti) 2 



dCi 



d(2J~ 

W2 VS. 1 



1 



1 



C2) 2 



In 



(wi - w 2 )(wi 



w 2 . 



47T 2 \(wi - W2)(wi - W2] 



£(wi, W 2 ) 



where Q is the Green function (1.5). With this step we have completed the proof 
of Theorem O 



5.4. Convergence to GFF: proof of Theorem 1.3 Our aim now is to prove 
the weak convergence of t/itHn (x, 77) (viewed as a generalized function on T>) to 
the w-pullback of the Gaussian free field GFF on H (see %\ 1.6- 1.7). In order to do 
that, we need an additional estimate: 



Lemma 5.6. For any e > and any s points (xi, T)i), • 
sarily pairwise distinct) we have the bound 

E(H N ( X i,Vi) ■ ■ - H N ( X s,Vs)) = 0(N*), 



, (Xs,Vs) € 2? (not 



N -> 00. 
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Proof. If all the points are distinct, we have a better bound O(l) by Theorem 
1.2 Next, assume that, say, (xi^Vi) = {X2, r l2)- Connect {x\,Vi) with the lower 
left edge by two paths which are close to each other only in a neighborhood of 
(xi, 771). As explained in <|3j we calculate (Hn(xi> Vi)) 2 as a product of sums over 
these two paths. Fix small S > 0, and in each sum consider separately the N^ +s 
terms corresponding to A~ s+^-neighborhood of (xitVi)- All other terms give a 
contribution of order O(l) because they involve points which are far apart, and so 
one can argue similarly to Proposition |4.11| and Theorem |1.2| On the other hand, 
the terms corresponding to close points are estimated using Proposition |4.10| and 



Lemma 



4.13 



The former gives gr owth of order O(lniV) 2 coming from the single 
integral in the correlation kernel (4.7), and the latter provides a bound 0(N 2S ) 
which comes from the double integral in (4.7). This concludes the proof. See also 
the end of Section 7 in [Ke n08| . □ 



Now we finish the proof of Theorem 1.3 We argue similarly to |BF08[ §5.5]. It 
suffices to establish that (see (1.9)) 



lim 7T 5/2 E TT 

N-yao VXXj=l 



MXi,Vt) H N(Xt,Vt)dXidm 



V 



(5.8) 



E 



FT , /0 4 (w- 1 (z 4 ))J(z t )GFF(z,)|dz l 



where cj>i, . . . ,<p s are smooth compactly supported test functions on T>. This conver- 
gence of moments implies the weak convergence ( 1.9 1 of y/ir J v 4>(x, v)Hn(Xi il)dxdr) 
to the corresponding Gaussian random variable J H </>(w _1 (z)) J(z) GFF(z)|rfz| 2 , as 
well as an obvious multidimensional analogue of this fact involving convergence to 
a Gaussian vector. 

is equal to 



The left-hand side of (|5. 



We split this integration over 
ciently far apart: 



(^(w- 1 ^)) J(zi)) E (iTj^w- 1 ^!)) . . . JMw" 1 ^))) . 

into two parts, one where the points Zj are suffi- 



H s , :={(zi,..., 

and the remaining part H s 
small and fixed. 

For the integration over 



z 5 ) G H 5 : \zi - Zj\ > N~i +S , 1 < i < j < s}, 
\ where some of them are close. As usual, 5 > is 



we use Propositions 5.3 and 5.4 



to write 



r 5 

/ TT \d Zl \ 2 (<j>i(™~\zi))J(zi)) E (Hn^- 1 ^)) . . . tf^w- 1 ^)) 



n\dz 



(w" 1 ^)) J(zi)) E ( GFF(zx) . . . GFF(z 5 )) + 0(N~l) 



Since the logarithms in E (GFF(zi) . . . GFF(z s )) (see f 1.6 1 are integrable around zero, 
we may replace the above integral over ET| by the same integral over 



The integral over the complement H 5 \ W s s is bounded using Lemma 5.6 



/ TT \d Zl \ 2 (^(w- 1 ^)) J(zi)) E (^(w- 1 ^)) . . . ^(w" 1 ^)) 
Jw\w s i=1 ^ 
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< const • (N-i +s ) 2 N*, 



where the constant depends only on our test functions <j)j . This last estimate implies 



the desired convergence (5.8), and completes the proof of Theorem 1.3 
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